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ROLE OF EJECTA CLUMPING AND BACK-REACTION OF ACCELERATED COSMIC RAYS IN THE 

EVOLUTION OF TYPE lA SUPERNOVA REMNANTS 

S. Orlando^ F. Bocchino^ M. Miceli^'\ O. Petruk^'*, M.L. Pumo'^''' 

ABSTRACT 

We investigate the role played by initial clumping of ejecta and by efficient acceleration of cosmic rays 
(CRs) in determining the density structure of the post-shock region of a Type la supernova remnant 
(SNR) through detailed 3D MHD modeling. Our model describes the expansion of a SNR through 
a magnetized interstellar medium (ISM), including the initial clumping of ejecta and the effects on 
shock dynamics due to back-reaction of accelerated CRs. The model predictions are compared to the 
observations of SN 1006. We found that the back-reaction of accelerated CRs alone cannot reproduce 
the observed separation between the forward shock (FS) and the contact discontinuity (CD) unless 
the energy losses through CR acceleration and escape are very large and independent of the obliquity 
angle. On the contrary, the clumping of ejecta can naturally reproduce the observed small separation 
and the occurrence of protrusions observed in SN 1006, even without the need of accelerated CRs. 
We conclude that FS-CD separation is a probe of the ejecta structure at the time of explosion rather 
than a probe of the efhciency of CR acceleration in young SNRs. 

Subject headings: cosmic rays — magnetohydrodynamics (MHD) — instabilities — shock waves — 
ISM: supernova remnants — supernovae: individual: SN 1006 



1. INTRODUCTION 

Nowadays it is widely accepted that supernova rem- 
nants (SNRs) are the site where cosmic ray (CR) dif- 
fusive shock acceleration occurs. Observations in var- 
ious bands support this picture through the detection 
of non-thermal emission that is compatible with being 
synchrotron or Inverse Compton radiation from CR elec- 
trons. Unfortunately, the direct evidence of CR ions in 
SNRs is difficult to find because they do not radiate effi- 
ciently. On the other hand, different indirect signatures 
of the presence of CR ions are largely discussed in the 
literature. The most popular is probably the separation 
between the forward shock and the contact discontinuity 
that has been me asured in young S NRs (e.g. SN1006 , 
Miceh et al.ll2009l and Tych o's SNR. iWarren et al.l[2005l : 
Cassam-Chenai et al.l 120071 ). In fact, current theories 
predict that a significant fraction of the energy of super- 
nova remnant shocks is channeled into CRs, determining 
modifications of the shock dynamics that depend on the 
efficiency of acceleration and injection processes of high 
energy particles. In particular, this energy losses would 
lead to a greater shock compression ratio and, as a conse- 
quence, to a thinner shell of shocked interstellar medium 
(ISM). 

An example of SNR in which the observed features 
have been interpreted as a consequence of the energy 
losses to CRs at the forward shock is SN 1006. In this 
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remnant, the observations have shown that the azimuthal 
profile of the ratio of the forward shock radius to the 
contact discontinuity radius Rfs/Rcd is fairly uniform 
(although very noisy) an d much lower than predicted 
for a non- modified sh ock (jMiceh et al.ll2009[ ). Recently 
iRakowski et al.l ()2011h have found and analyzed clumps 
of ejecta close to or protruding beyond the main blast 
wave of SN 1006 that have been interpreted in the con- 
text of an upstream medium modified by the saturated 
nonresonant Bell instability which enhances the growth 
of Rayleigh- Taylor (RT) instabilities at the contact dis- 
continuity. 

However, some pieces of evidence are now accumulat- 
ing that are difficult to explain in ter ms of acceleration 
of CR particles. Some authors (e.g. iBlondin fc EUisonl 
1200 ll : rWandl2011l and references therein) noted that ex- 
treme energy losses to accelerate the CRs are needed to 
allow a significant fraction of the ejecta to approach or 
even overtake the forw ard shock, th us explaining the thin 
shell of shocked ISM. IWand (|2011t ) analyzed the evolu- 
tion of RT instabilities in Type la SNRs undergoing CR 
particle acceleration and found that, even with very effi- 
cient acceleration of CRs (i.e. assuming an effective adia- 
batic index 7cff ~ 1-1), significantly enhanced mixing and 
perturbation of the remnant outl ine are not expected. A 
similar conclusion was reached bv lFraschetti et all (|2010[ ) 
who found that the development of RT instabilities in 
SNRs is not drastically affected by CR particle acceler- 
ation. In addition, these studies suggest that the high 
occurrence of protrusions in young SNRs is not the con- 
sequence ofRTjnstaMitifis by accelerated CRs 
(see also lWang fc Chevalieill2001[ ). Another evidence dif- 
ficult to explain in terms of acceleration of CR particles is 
the ratio Rfg/Rcd measured in SN 1006 that is lower than 
predicted by non-modified shock models even in regions 
dominated by thermal emission where the CR accelera- 
tion efficiency is supposed to be low (e.g. iMiceli et al.l 
I2009D . All these studies cast some doubts on whether 
the back-reaction of accelerated CRs is the main respon- 
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sible of the enhanced intershock instabihties observed in 
young SNRs (e.g. SN 1006 and Tycho's SNR). 

On the other hand, spectropolarimetric studies of 
SNe la show the presence of asymmetries with different 
magnitude and orientation for different elements in the 
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All these features have been interpreted as being due 
to c lumpy structures in the outer layers of the ejecta 
(see iHole et all |2010| and references therein) and some 
authors suggested that ejecta clumps of intermediate- 
mass elements ca,n be f orged in the explosi on of SN la 
(e.g. IWang et al.l 120031: ILeonard et al.l 120051 ) or may be 
due to the interaction of the ejecta with a dense, c lumpy, 
and d isk-like circ umstellar envi r onmen t (e.g. Wang et al.l 
I2004D . Recentlv [Maeda et al.l (|2010aD have shown that 
asymmetries in t he explosion can be a generic feature in 
SNe la (see also iMaeda et al.l 120109 ) . and these asym- 
metries, in turn, may lead to a clumpy structure of the 
ejecta. In the light of these considerations, it is there- 
fore important to investigate the role of ejecta clumping 
on the evolution and morphology of Type la SNRs. In 
particular we wonder whether the thermal and density 
structure of the post-shock region of a young SNR orig- 
inates mainly from the clumpy structure of the ejecta 
rather than as a consequence of back-reaction of accel- 
erated CRs. The density inhomogeneities in the ejecta 
can enhance the growth of RT instabilities, causing the 
ejecta material to move closer to the main blast. The 
question is: can the ejecta clumping enhance the growth 
of RT instabilities up to a level that allows clumps of 
ejecta to reach and possibly overtake the forward shock? 

Here we investigate this issue by developing a three- 
dimensional (3D) MHD model describing the expansion 
of a SNR through a magnetized medium, including, for 
the first time, the (non-uniform) ambient magnetic field, 
the initial ejecta clumping, and the effects on shock dy- 
namics due to back-reaction of accelerated CRs. The 
paper is organized as follows: in Sect. [5] we describe the 
MHD model and the numerical setup; in Sect. |3|we de- 
scribe the results and, finally, we draw our conclusions 
in Sect. [H 

2. MHD MODEL AND NUMERICAL SETUP 

The evolution of a SNR can be characterized by dis- 
tinct stages depending on the physical pr ocess domi- 
nating its dynamics (e.g. iChevalieil I1977D . This pa- 
per focuses on young SNRs, i.e. remnants that have 
evolved from the ejecta-dominated stage through the 
Sedov- Taylor stage. Pioneering comprehensive studies of 
the dynamics of these remnants, preceding the onset of 
dynamically significant radiative losses and/or pressure 
confinement b y the ambient medium, are given in the 
literature (e.g. iMansfield fc Sal peteilll974l : iFranco et al.l 
119941: iTruelove fc McKcc 1999) and are mostly based on 
analytic and numerical ID hydrodynamic models. Sub- 
sequently several 2D and 3D hydrodynamic and MHD 
models describing the evolution of the remnant through 
the ISM have been developed. 

Here we adopted the 3 D MHD model discussed by 
lOrlando et al.l (j2007ll201lD . extended to describe the ini- 
tial ejecta clumping and to include the effect of larger 
compressibility of plasma around the shock due to the 



back reaction of accelerated CRs. The shock propagation 
is modeled by numerically solving the time-dependent 
ideal MHD equations of mass, momentum, and energy 
con servation in a 3p car tesian coordinate system (x, y, z) 
fsee lOrlando et aLll2007l for details). In order to trace the 
motion of the ejecta material and study its dynamics, we 
considered a passive tracer associated with the ejecta. 
The continuity equation of the tracer is solved in addi- 
tion to our set of MHD equations; the ejecta material is 
initialized with Coj = 1, while Ccj = in the ISM. The 
calcu lations were performed using flash (jFrvxell et al.l 
l200Clf ). an advanced multi-dimensional MHD code for as- 
trophysical plasmas, including the adap tive mesh refine- 
ment through the PARAMESH library ()MacNeice et al.l 
l2000f ). and extended with additional computational mod- 
ules to handle the back-reaction of accelerated CRs. 

The effects of shock modification are i ncluded in the 
MHD model by following the approach of iFerrand et al.l 
(|2010f) and extending their method to MHD models. In 
particular, our model includes an effective adiabatic in- 
dex 7off which depends on the injection rate rj of parti- 
cles (i.e. the fraction of ISM particles entering the shock 
front). The adiabatic inde x on the shock is var ied due to 
particle acceleration as in lElhson et al.l (l200l ( see also 
iFerrand et al.|[2010l ). At each time-step of integration, 
the adiabatic index is calculated at the shock front and 
then is advected within the remnant, rema ining constant 
in eac h fluid element. As discussed by IFerrand et al.l 
(|2010( l. the latter assumption implies that each fluid ele- 
ment remembers the effect of shock modification induced 
by particle acceleration at the time it was shocked. 

For the purposes of the present paper, we assume 
that the maximum injection rate rj is large enough (e.g. 
rj w 10^"^, namely when shock modifications are strong 
and immediate) so that the effective adiabatic index 
at the initial conditions of our simulations has already 
reached its minimurn value and slightly depends on time 
(jFerrand et al.ll2010D . We assume therefore the effective 
adiabatic index not depending on time and consider its 
minimum value 7inin as a free parameter. On the other 
hand, the injection rate is expected to depend on the 
shock obliquity (i.e. the angle between the unperturbed 
external magnet ic field and the normal to the shock; e.g. 
iVolk et al.|[2003D . We allow therefore that the effects of 
shock modification on the fluid dynamics (and, there- 
fore, the effective adiabatic index) vary in space as a 
function of the obliquity angle. We assume no magnetic 
field amplification due to CRs, and no back-reaction of 
accelerated CRs at the reverse shock, although the sug- 
gestion that CR particles can be efficiently accelerated 
also at th e reverse shock is la rgely debated in the litera- 
ture (e.g. lEUison et al1 l2005^ 

The index 7eff is calculated at the shock front by us- 
ing a parametrized function depending on the obliquity 
angle Oo and characterized by a parameter representing 
the minimum value of the adiabatic index 7niin that is 
possible to reach during the simulation: 



7cff = 7 - (7 - 7min) X A(eo) 



(1) 



where 7 = 5/3 is the adiabatic index and /<;(0o) is 
a function defined in the range [0, 1] depending on the 
obliquity angle Oo and describing the variations of 7off 
over the surface of the remnant shock. In analogy with 
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the description given by iFulbright fc Reynolds! (|1990[ ) 
for the quasi-parallel, quasi- perpendicular, and i s otropi c 
injection models (see also I Orlando et alj l2007l |2Q11D . 
we model the variations of 7eff over the shock surface 
through the functions /<;(6o) = cos^ Og (i.e. jcB is min- 
imum at parallel shocks), /?(6o) = sin 9s (7off is min- 
imum at perpendicular shocks), and A(6o) = 1 (7cff is 
uniform at the shock front and equal to its minimum 
value), where Qg is the angle between the shock nor- 
mal and the post-shock magnetic field and is related 
to 6o by the expression cosOg — cr~^ cos 0o and a 
is the shock compression ratio. The first case follows 
the quasi-parallel injection scenario, leading to a 3D 
polar-caps structure of the remnant, whereas the second 
and the third cases follow the quasi-perpendicular and 
isotropic injection models, respectively, producing a 3D 
equatorial-belt structure of the remnant. Note that the 
third case (isotropic) is intended to be the extreme case 
in which the shock modification is the largest everywhere 
at the forward shock with no obliquity dependence. 

As for the density structure of the ejecta, we inves- 
tigated: the exponential profile that has been shown 
to be the most represent ative of explosion mode ls for 
thermonuclear SNe (Dwar feidas fc Chevalieij|l998f) . and 
the power-law profile with index n = 7 that has been 
used to represent d efiagration models (jChevalieil Il983t 
iNomoto et aH 119841 ) ■ We also assume that the initial 
ejecta has a clumpy structure. The clumps have been 
modelled as per-cell random density perturbations^ de- 
rived from a power-law probability distribution^ (index 
n — —1) that is characterized by a parameter z^max rep- 
resenting the maximum density perturbation allowed in 
the simulation. Figure [T] shows the power-law probabil- 
ity distributions of the perturbations used in this paper 
for the two ejecta density profiles considered. We ex- 
plored maximum density perturbations ranging between 
1.5 and 5; we explored density clumps of ejecta with size 
either 1% or 2% of the initial diameter of the remnant 
DgnrO- As discussed in Sect. 13.21 initial clump size in the 
range explored here leads, after 1000 yr of evolution, to 
density features with characteristic size comparable to 
those observed in SN 1006. 

It is interesting to note that the range of clump size 
investigate d in this pa. p er is a lso in agreement with that 
derived by iHole et al.l ()2010[ ) for SNe la. In particular 
these authors compared the results of their semi-analytic 
code for modeling polarized line radiative transfer within 
3D inhomogeneous rapidly expanding atmospheres with 
spectropolarimetric observations; they found that the 
model reproduces the observed range of values of peak 
line polarization if the clumps have radius in the range^ 

The density perturbation of each clump is calculated as the 
ratio of the mass density of the resulting clump to the local average 
density in the region occupied by the clump if the perturbation was 
not present. 

* Since no observational clues are available on the distribution of 
density perturbations of the clumps, it is reasonable to assume that 
most of the clumps are characterized by small density perturbations 
and few of them by large perturbations. To this end, for the sake 
of simplicity, we assume a power-law probability distribution with 
index n = — 1. 

^ Sizes are given in velocities because of the linear dependence 
of velocity on distance in the homologous flow of SN ejecta (this 
structure is som etimes referred to as a pseudo-Hubble flow; see 
IHole et al.ll2010h . 
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Fig. 1. — Probability distribution functions of the random per- 
turbation of mass density of the clumps for the two ejecta density 
profiles considered in this paper: the exponential profile (upper 
panel) and the power-law profile with index n = 7 (lower panel). 
The density perturbation of each clump is calculated as the ratio 
of the density of the resulting clump to the local average density 
in the region occupied by the clump if the perturbation was not 
present. 

1000 — 6000 km s^^. At the time of our initial condition 
(« 10 yr since the SN explosion) , the e ffective range of 
clump size derived by iHole et all ()2010l ) corresponds to 
0.016—0.13 pc to be compared with the size of the clumps 
modelled here ranging between 0.01 and 0.02 pc. As an 
example. Fig. [5] shows the initial spatial distribution of 
ejecta clumps for a model with the highest density per- 
turbation and largest clump size. A summary of all the 
simulations discussed in this paper is given in Table [T] 

Note that the ejecta clumps are presumably relics of 
the deflagration of the outer layers of the exploding 
star (as suggested by theoretical arguments and observa- 
tions). In principle, therefore, the clumps are expected 
to be concentrated in a shell within the ballistically ex- 
panding ejecta rather than being distributed in the whole 
unshocked ejecta as done here. On the other hand, in our 
simulations, the ramp profile of the initial velocity of the 
ejecta makes the clumps in the outer layers those with 
the highest speed, so that the shocked ISM is mostly 
perturbed by such clumps. Concerning the focus of this 
paper, namely the structure of the RT mixing in the re- 
gion between the forward and reverse shocks, we do not 
expect therefore significant changes to the our results if 
considering a distribution of clumps concentrated in the 
outer layers of the ejecta. 

As initial conditions, we adopted parameters appro- 
priate to reproducing the SNR SN1006 after 1000 yr 
of evolution: we assumed an initial spherical remnant 
with radius i?fso = 0.5 pc (corresponding to an initial 
age of ~ 10 yr), originating in a progenitor star with 
mass of 1.4 Mgun, and propagating through an unper- 



4 S. Orlando et al. 



TABLE 1 

Adopted parameters and initial conditions for the MHD models of the SNR 



Model 


eject a 


ejecta 


shock 


a 

t^max 


clump 


injection 


min(7cfl) 


initial 


abbreviation 


profile 


clump. 


mod. 




sizc^ 


efficiency 




age [yr] 




EXP*^ 


no 


no 








5/3 


10 


REF-PL 


FLAW 


no 


no 








5/3 


10 


h/X-Ol.o-Dl 


EXP 


yes 


no 


1.5 


l7o 




5/3 


10 


TT>"V' Ol c "Hin 

ii/X-O1.0-U2 


EXP 


yes 


no 


1.5 


00/ 




5/3 


10 


il/A-Oz.5-Dl 


EXP 


yes 


no 


2.5 






5/3 


10 


il/A-Cz.o-JJz 


T^VT) 


yes 


no 


2.5 






0/3 


10 


il/A-Co.U-JJi 


T^VT) 

il/AF 


yes 


no 


5.0 


1 0/ 
i /o 




0/3 


10 


ii/A-Oo.U-U2 


EXP 


yes 


no 


5.0 


00/ 




5/3 


10 


FL-Cl.o-Dl 


PL AW 


yes 


no 


1.5 


1 0/ 

l7o 




5/3 


10 


FL-Cl.o-Di 


PL AW 


yes 


no 


1.5 


oO/ 




5/3 


10 


rLi-K^Z.O-U L 


T>T AIS^ 

tLAW 


yes 


no 


2.5 


1 0/ 
i /o 




0/3 


10 


T)T /^O C T\0 


PL AW 


yes 


no 


2.5 


00/ 

i7o 




5/3 


10 




T>T ATXJ 

tLAW 


yes 


no 


5.0 


1 0/ 

i /o 




0/3 


10 


T)T /^c n T^O 


PL AW 


yes 


no 


5.0 


00/ 




5/3 


10 


A- Q F A H- G 1 . 1 


EXP 


no 


yes 






QPAR*^ 


1.1 


10 


TTV OPAR r^l Q 

11; A- T rt. IX- O i • O 




no 


yes 








4/3 




PL-QPAR-Gl.l 


FLAW 


no 


yes 






QPAR 


1.1 


10 


PL-QPAR-G1.3 


FLAW 


no 


yes 






QPAR 


4/3 


10 


EX-ISO-Gl.l 


EXP 


no 


yes 






ISO-'' 


1.1 


10 


PL-ISO-Gl.l 


FLAW 


no 


yes 






ISO 


1.1 


10 


EX-C3.5-D1-QPAR-G1.3 


EXP 


yes 


yes 


3.5 


1% 


QPAR 


4/3 


10 


EX-C3.5-D2-QPAR-G1.3 


EXP 


yes 


yes 


3.5 


2% 


QPAR 


4/3 


10 


PL-C3.5-D1-QPAR-G1.3 


FLAW 


yes 


yes 


3.5 


1% 


QPAR 


4/3 


10 


PL-C3.5-D2-QPAR-G1.3 


FLAW 


yes 


yes 


3.5 


2% 


QPAR 


4/3 


10 


EX-C1.5-D1-2YR 


EXP 


yes 


no 


1.5 


1% 




5/3 


5 


EX-C5.0-D2-2YR 


EXP 


yes 


no 


5.0 


2% 




5/3 


5 



" Maximum perturbation of mass density; percentage of the initial diameter of the remnant; exponential profile; power-law profile; 
quasi-parallel injection; f isotropic injection. 

The initial total energy Eq = 1.5 x lO'^^ ergs leads to 
a remnant radius Rsnr ~ 8.5 pc at < = 1000 yr and is 
partitioned so that > 99% of the SN energy is kinetic. 
The remnant expands through a homogeneous isothermal 
medium of plasma number density n = 0.05 cm~'^ and 
temperature T — 10^ K. The initial a mbient magnetic 
field configuration is that suggested by iBocchino et alj 
(|2011l) for SN 1006 and resulting from the comparison 
of radio observations of SN 1006 with MHD models: the 
ambient magnetic field is characterized by a non-zero gra- 
dient of its strength perpendicular to the average mag- 
netic field that leads to a variation of \B\ of about a 
factor 1.4 over a scale of 10 pc. In all our simulations, 
the magnetic field strength is ~ 3 fiG in the environment 
of the explosion site. We follow the remnant evolution 
for 1000 yr. 

The computational domain extends 24 pc in the x, y, 
and z directions. Special emphasis was placed on captur- 
ing the enormous range in spatial scales in the remnant. 
To this end, we exploited the adaptive mesh capabilities 
of the FLASH code by using 11 nested levels of resolution, 
with resolution increasing twice at each refineme nt level. 
The refinement /derefinement criterion adopted (|Lohneil 
|1987[ ) follows the changes in mass density, temperature, 
and tracer of ejecta. In addition, the calculations were 
performed using also an automatic mesh derefinement 
scheme in the whole spatial domain that kept the com- 
putational cost approximately constant as the blast ex- 
panded: the maximum number of refinement levels used 
in the calculation gradually decreased from 11 (initially) 
to 7 (at the final time) following the expansion of the 
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Fig. 2. — Left panel: initial spatial distribution of plasma density 
along the x-axis for a model either with (red line; run EX-C5.0- 
D2 in Table [TJ or without (black line; run REF-EXF) the ejecta 
clumping. In both models the total mass of ejecta (integrated 
over the whole volume) is 1.4 Msun (see text). Right panels: initial 
spatial distributions of ejecta clumps with density perturbation in 
the range either [3.5 — 4] (upper panel) or [4.5 — 5] (lower panel) 
in run EX-C5.0-D2. 

turbed magneto-static medium. Note that we payed 
particular attention to have an initial total mass of ejecta 
Mej = 1.4 Msun in all the simulations considered here, 
including those with a clumpy structure of the ejecta. 
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blast and keeping roughly the same number of grid zones 
per radius of the remnant. At the beginning (at the end) 
of the simulation, this grid configuration yielded an ef- 
fective resolution of « 2.9 x 10"'^ pc (« 4.6 x 10~^ pc) at 
the finest level, corresponding to w 170 zones per initial 
radius of the remnant (~ 190 zones per final radius of 
the remnant). The effective mesh size varied from 8192'^ 
initially to 512^ at the final time. 

We also performed two additional simulations with the 
same parameters of runs EX-C1.5-D1 and EX-C5.0-D2 
but starting as early as « 2 yr after the SN explosion (the 
initial spherical remnant has radius i?fso = 0.125 pc) to 
check if the results depend on the time when the clumpy 
structure of the ejecta is initialized. The results of this 
comparison are discussed in Appendix [XI 

3. RESULTS 

3.1. Effects of back-reaction of accelerated cosmic rays 

As a first step, we analyzed the effects of back-reaction 
of accelerated CRs on the separation between the blast 
wave and the contact discontinuity, by considering mod- 
els accounting for the shock modification by accelerated 
CR particles but without initial clumping of ejecta^°. A 
recent comprehensive study of these effects on the de- 
velopment of RT instabilities in young SNRs is given by 
IWand ([20111 ) (see also references therein). Our study 
differs from previous works in that it includes magnetic 
fields and a possible dependence of the CR particle ac- 
celeration on the obliquity angle. In particular, we fo- 
cused on the isotropic and quasi-parallel scenario dis- 
cussed in Sect. [21 the results for models assuming quasi- 
perpendicular injection are expected to be analogous to 
those discussed here for quasi-parallel injection, showing 
a modulation of the shock modification with the obliq- 
uity angle. 

As expected for cases in which the magnetic field has a 
comp onent parallel to the s urface of the contact disconti- 
nuity ([Chandrasekha J 196 1[ ) , the magnetic field limits the 
growth of hydrodynamic instabilities through the tension 
of field lines which maintain a more laminar flow around 
the contact discontinuity. The energy losses to CRs at 
the forward shock lead to a greater sho ck compression ra- 
tio in all the cases examined (see also [Blondin fc EUisonl 
[200lt [mMl[20Tll) . As a consequence, the density of 
the shocked ISM is greater and the separation between 
the blast wave and the contact discontinuity is shorter 
than predicted for a non-modified shock in regions with 
7efi < 5/3, i.e. where the back-reaction of accelerated 
CRs is efficient. In the quasi-parallel case, since the 
back-reaction of CRs is more effective at parallel shocks, 
the shock modification is modulated with the obliquity 
angle. As an example of this case. Fig. [3l presents the 
results for a model with an exponential profile of the 
initial ejecta density after 1000 yr of evolution (run EX- 
QPAR-Gl.l; see Table [TJ. In this model we also assumed 
extreme energy losses to accelerate the CRs, so that the 
minimum effective adiabatic index is j^ff — l-l- The 
modulation of the back-reaction of accelerated CRs with 
the obliquity angle is evident in the figure, showing a 
larger compressibility and higher values of plasma den- 

Note that, in these simulations, we did not introduce any seed 
perturbation. The departures from spherical symmetry are entirely 
due to the mesh and to possible numerical fluctuations. 
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Fig. 3. — 2D sections in the {x,z) plane of the spatial distribu- 
tion of plasma number density (upper panel) and effective adiabatic 
index (middle panel), after 1000 yr of evolution, for a model ac- 
counting for the shock modification by accelerated CRs but without 
initial clumping of ejecta (run EX-QPAR-Gl.l). The lower panel 
shows the 3D volume rendering describing the spatial distribution 
of the effective adiabatic index. The index ■y^ff is minimum in red 
regions (sec colour bar). The white lines in the upper two pan- 
els and the yellow lines in the lower panel are sampled magnetic 
field lines. The violet surface in the lower panel tracks the ejecta 
material. 

sity at parallel shocks. Such a modulation is absent in 
the isotropic case where the effects of CR particle ac- 
celeration are the same everywhere at the shock front 
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(runs EX-ISO-Gl.l and PL-ISO-Gl.l). In these cases 
the plasma compressibility is the largest everywhere at 
the shock front, and the post-shock magnetic field can 
reach values up to « 50 — 70 /iG at perpendicular shocks. 
It is worth mentioning that in both the quasi-parallel and 
isotropic cases, the simulations do not show any signifi- 
cant perturbation of the remnant outline and occurrence 
of protrusions after 1000 yr of evolution, even assum- 
ing extreme energy losses to accelerate the CRs. These 
results are in agreement with previous studies showing 
that enhanced RT mixing due to efficient particle ac- 
celeration determines only a slight perturbation of the 
forward shock near the epoch of young SNRs as SN 1006 
or Tycho (e.g. iBlondin fc Ellisonll200ll: IWanell201ll ). 

We investigated the effect of accelerated CRs on the 
separation between the blast wave and the contact dis- 
continuity, by deriving the azimuthal profiles of the ratio 
of the forward shock radius to the contact discontinu- 
ity radius Rfs/Rcd from the models. The position of the 
forward shock was estimated from 2D maps of projected 
emission-measure- weighted temperature (T) as the jump 
in (T) in the direction of compression (determined by 
looking at the velocity field) at temperatures T > 1 MK. 
The position of the contact discontinuity was estimated 
by using the passive tracer Coj included in the model 
(see Sect. [2]): during the remnant evolution, the ejecta 
and the shocked ISM mix together, leading to regions 
with < Ccj < 1; at any time t the density of ejecta 
material in a fluid cell is given by poj = pCcy We de- 
rived the position of the contact discontinuity from 2D 
maps of projected poj as the local peak of pcj closest to 
the forward shock in the direction of compression. The 
azimuthal profiles derived from the models in such a way 
are directly comparable with obse rvations and, i n par tic- 
ular, with the profiles derived bv iMiceli et al.l (|2009[ ) in 
the ana lysis of the observations of SN 1006 (see Fig. 6 in 
IMiceli et all (20091 ) ■ 

Fig. m shows the azimuthal profiles of Rfs/Rcd derived 
from the models when the aspect angle is 90° (i.e when 
the average magnetic field is perpendicular to the line of 
sight). The green line is the profile derive d from the ob- 
servations of SN1006 (jMiceli et al.ll2009t) and the black 
lines are the reference models with no shock modifica- 
tion and no ejecta clumping. The gray stripes mark the 
regions where, in the models, the acceleration of CRs 
is the largest. In the quasi-parallel case, we found that 
the modeled profiles are modulated by the obliquity an- 
gle and, in general, are higher than those observed (see 
red and blue lines in Fig. \^ . The observations can be 
reproduced only in limited regions where the effect of ac- 
celerated CRs is the highest. Aspect angles lower than 
90'^ make the comparison between models and observa- 
tions worse because the regions of efficient CRs accelera- 
tion would not be at the limb (so that the ratio Rfs/Rcd 
increases). On the other hand, the models with no obliq- 
uity dependence of the back-reaction of accelerated CRs 
and 7cff = 1.1 reproduce the observed profiles quite well 
(magenta lines in Fig. |4]). These results suggest that 
the observations could be reproduced only if the back- 
reaction of accelerated CRs is extreme (i.e. jeS ~ 1.1) 
and independent of obliquity angle (i.e. the CRs accel- 
eration and escape should be ubiquitous at the forward 
shock). 
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Fig. 4. — Azimuthal profiles of the ratio of tlie forward sfiock ra- 
dius to the contact discontinuity radius Rfs/ R^d for models without 
ejecta clumping and with initial ejecta density profile either expo- 
nential (upper panel) or power-law (lower panel). The black lines 
mark the profiles derived from the reference models not accounting 
for the back-reaction of accelerated CRs (runs REF-EX and REF- 
PL). The red and blue lines mark the profiles derived from models 
including the shock modification modulated by the obliquity angle 
either with minimum ■y^fi = 4/3 (red lines; runs EX-QPAR-G1.3 
and PL-QPAR-G1.3) or minimum ■y^g = 1.1 (blue lines; runs EX- 
QPAR-Gl.l and PL-QPAR-Gl.l). The magenta fines mark the 
profiles derived from models including the shock modification with 
no obliquity dependence and ■y^g = 1.1 (runs EX-ISO-Gl.l and 
PL-ISO-Gl.l). The greenline marks the profile derived from the 
observations of SN 1006 jMiceli et al. 2009). The gray stripes mark 
the regions where the acceleration of CRs is the largest. 



3.2. Effects of ejecta clumping and instability 

As a next step, we investigated the effects of ejecta 
clumping on the evolution and morphology of the rem- 
nant by considering models without back-reaction of ac- 
celerated CRs and accounting only for the ejecta clump- 
ing. In addition to the spectrop olarimetric studie s of 
SNe discussed in Sect. [1] (see also iHole et all 120101 and 
references therein), a widespread dumpiness of ejecta is 
also suggested by X-ray and radio observations, showing 
knots located near the edge of the remnants, and out- 
ward protrusions in many cases surrounding the knots 
(e.g. iHwang fc Gotthell 119971: IVelazguez eFal] 119981: 
iRakowskiet al.l 1201 iD . All these features cannot be ex- 
plained by instabilities generated by linear perturba- 
tions and have been interpreted as being due to clumps 
of ejecta expanding in to the intershock region (e.g. 
I Wang fc Chevalieii[200ll ). The interactions among the 
clumps of ejecta are expected to contribute to seed the 
RT instabilities and enhance their growth, thus strongly 
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influencing the final morphology of the remnant. 

The basic physics of the evolution of a single clump 
of ejecta expanding through the intershock structure of 
a SNR is similar to that for the interaction o f a shock 
with a cloud of the ISM (e.g. 'Klein et al. 1994) and has 
been extensively discussed bv iWang fc Chevalier (200 li ). 
The major factors in the clump-remnant interaction are 
the density contrast of the clump with respect to the 
ISM, the clump size, and the position of the clump in 
the initial distribution of ejecta (or, alternatively, the 
time of initiation of the clump-shock interaction). In 
general, after passing through the reverse shock, the sin- 
gle clump evolves toward a core-plume structure with 
a crescent-like shape characterized by Kelvin-Helmholtz 
(KH) instabilities developing in the downstream region. 
As the clump travels through the intershock structure, 
RT instabilities develops on the upstream side of the 
clump, leading to its progressive fragmentation. Depend- 
ing on its initial density contrast, size, and time of ini- 
tiation of the clump-shock interaction, the clump can 
reach the forward shock, causing a bulge on the remnant 
outline as the ram pres sure pushes material ahead (see 
IWang fc Chevaliei1l2001l for a detailed description). This 
is the way ejecta protrusions form. After the clump is 
completely fragmented, the bulge (the protrusion) dis- 
appears and the clump fragments are mixed with the 
shocked ISM and swept back in the remnant. The per- 
turbation of the forward shock front by the interaction 
with the clumps is more likely during the early phases of 
the remnant evolution when the density contrast between 
the ejecta clumps and the ISM is larger. 

In our case, we are assuming that the ejecta structure is 
formed by hundreds of thousands of clumps modelled as 
per-cell random density perturbations (see Sect. [5] and 
Fig. [5]); in each simulation, the clumps have the same 
size and are characterized by different density contrasts 
(i.e. different density perturbations) and different posi- 
tions (i.e. the time of initiation of the interaction of each 
clump with the reverse shock is different). The clump- 
remnant interaction therefore is complicated by the mul- 
tiple interactions among clumps with different density 
contrast and velocity. In addition, our model includes 
the magnetic field which is known to limit the growth 
of hydrody namic instabilit i es in t he shock-cloud int erac- 
tion (e.g. iMac Low et"all 119941: Ijones et al.lll996[) due 
to the tension of the magnetic field lines which maintain 
a more laminar flow around the cloud surface (see also 
iFragile et ani2005HOri"ando et al.ll2008| ). In the present 
case, during the clump evolution, the magnetic field is 
expected to be trapped at the nose of the clump, leading 
to a continuous increase of the magnetic pressure and 
field tension there that limit the growth of RT instabili- 
ties responsible for the clump fragmentation. As a result, 
the clumps are e xpected to survive for a lo nger time than 
those studied bv lWang fc Chevalieii (|2001[ ) (their simula- 
tions do not include the magnetic fleld), increasing their 
probability to reach the forward shock. 

As an example. Fig. [5] shows a close-up view of the 
remnant limb for the model EX-C5.0-D2, illustrating the 
magnetic field strength (upper panel) and the plasma 
density distribution (lower panel) aX t — 1000 yr. As 
expected, the magnetic field follows the plasma struc- 
tures formed during the evolution of the clumps with 
preferentially radial components around the RT fingers. 
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Fig. 5. — Close-up view of the remnant limb for the model EX- 
C5.0-D2 at t = 1000 yr, showing a colour-coded cross-section image 
of the magnetic field strength (G; upper panel) and a composite 
cross-section image (lower panel) combining the plasma density 
(cm~'^) of the shocked ISM (red) and that of the ejecta (blue- 
green). The contours enclose the cells consisting of the original 
ejecta material by more than 10%. The magnetic field is described 
by the superimposed arrows the length of which is proportional to 
the magnitude of the field vector. 

The magnetic field is strongly modified by the clumps 
and it can be enhanced by up to two orders of mag- 
nitude (\B\ K, 100 /iG, whereas the unperturbed mag- 
netic field strength is « 2.5 /^G) in the ejecta clumps 
(see red regions in the upper panel of Fig. [5]). Note 
that, in model EX-G5.0-D2, no back-reaction of acceler- 
ated CRs is taken into account and the magnetic field 
in inter-clumps regions at the forward shock {\B\ k, 
10 /iG) is that predicted for non-modified shocks, namely 
much lower than that measured in the X-ray rims of 
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Fig. 6. — 3D rendering of mass density for models with an exponential (upper panels) or a power-law (lower panels) initial profile of 
ejeeta and no shock modifieation by CRs, after 1000 yr of evolution. Left panels show the maps derived from the reference models not 
including the ejecta clumping (runs REF-EX and REF-PL). The initial clumpy structure of the ejeeta is characterized by clumps either 
with size fa 1% of the initial diameter of the remnant DsnrO and maximum density perturbation Vmax = 1.5 (center panels) or with size 
~ 2% of the initial diameter of the remnant and maximum density perturbation i^max = 5.0 (right panels). The white contours enclose the 
ejecta material (i.e. cells consisting of the original ejecta material by more than 50%). 

sity perturbation (i.e. going from the left to the right 
panel of Fig. both the perturbation of the remnant 
outline and the occurrence of ejecta protrusions increase. 
Fig. [6] also shows that the characteristic size of the den- 
sity features formed within the remnant is comparable to 
that of the features observed, for instance, in SN 1006. 

A remarkable feature of the simulations including the 
ejecta clumping is the occurrence of several protrusions 
due to clumps of ejecta overtaking the forward shock. 
Fig. [7] shows composite images of the SNR combining 
the square of plasma density of the shocked ISM (red) 
and that of the ejecta (green and yellow), both projected 
along the line-of-sight, for models EX-C5.0-D2 and PL- 
CS. 0-D2. The protrusions are evident in both cases and 
are due to clumps with high density contrast originating 
from the outer layers of the ejecta. Our calculations show 
that the number of protrusions aX t — 1000 yr is higher 
for larger size of the clumps and higher density contrasts 
of the clumps and decreases with the age of the remnant. 
In fact the simulations showed that, during the remnant 
evolution, new protrusions are continuously formed and, 
subsequently, disappear when the clumps responsible for 
them are decelerated and the forward shock front catches 



SN1006 (50 < 1^ 1 <150/iG;IBerezhko"eFan[2003L [200l 
lAcero et alllfoTot IPetruk et al.l[20Tll[20T2h . 

As examples. Fig. [S] shows the 3D rendering of plasma 
density for the reference cases without clumping (runs 
REF-EX and REF-PL) and for the hmit cases with 
clumping considered in this paper, namely models with 
an ejecta structure characterized either by clumps with 
small size and low density perturbations (runs EX-C1.5- 
Dl and PL-C1.5-D1 in Tabled]) or by clumps with large 
size and high density perturbations (runs EX-C5.0-D2, 
and PL-C5.0-D2). The figure shows that the enhanced 
intershock RT mixing can easily spread the ejecta mate- 
rial close to, or even beyond, the average radius of the for- 
ward shock, depending on the size and density contrast 
of the initial clumps. This can occur very soon after the 
explosion, depending again on the size and density con- 
trast of the clumps seeding the instabilities. As a result, 
we found that: 1) the RT mixing reaches the forward 
shock front possibly perturbing the remnant outline, 2) 
knots and filamentary structures characterize the rem- 
nant morphology, and 3) clumps of ejecta can be very 
close to or even protruding beyond the main blast wave 
leading to evident knots near the remnant edge as ob- 
served, for instance, in SN 1006 (Rakow ski et al. 2011) 
and Tycho's SNR (e.g. IVelazguezeTaLl Il998l) . In gen- 
eral, increasing the initial size of the clumps or their den- 



up with them (see also Wan g fc Chev alier 2002). In this 
process, the clumps contribute in the perturbation of the 
remnant outline and in the formation of plasma features 
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Fig. 7. — Composite images of the remnant limb showing the 
square of plasma density of the shocked ISM (red) and that of the 
ejecta (green and yellow), both projected along the line-of-sight, for 
models EX-C5.0-D2 (upper panel) and PL-C5.0-D2 (lower panel) 
after 1000 yr of evolution. 

in the outer part of the remnant. 

Also in this case, we compared the azimuthal profiles 
of the ratio Rfs/ Red derived from the models with that 
observed in SN 1006 (see Fig. H]). We found that the 
initial clumping of ejecta makes the azimuthal profiles of 
Rfs /Red fairly uniform and lower than expected for mod- 
els without a clumpy structure of the ejecta and compa- 
rable with models accounting for extreme and ubiqui- 
tous acceleration of CR particles at the forward shock 
(i.e. isotropic models with 7off ~ 1.1; compare Fig. 2] 
and Fig. In particular, we found that, in the case of 
SN 1006, the observed profile can be reproduced by mod- 
els with a maximum density perturbation of ejecta fmax 
ranging between 2.5 and 5, and with initial size of ejecta 
clumps of the order of 2% of the initial diameter of the 
remnant (see the blue lines in middle and lower panels 
in Fig.©. 

3.3. Ejecta clumping and cosmic rays acceleration 

As a last step, we have investigated the effects of back- 
reaction of accelerated CRs on the remnant morphol- 
ogy in the presence of ejecta clumping through simula- 
tions including both physical processes (runs EX-C3.5- 
D1-QPAR-G1.3, EX-C3.5-D2-QPAR-G1.3, PL-C3.5-D1- 
QPAR-G1.3, and PL-C3.5-D2-QPAR-G1.3 in Tabled]). 
We found that when the CR acceleration efficiency de- 
pends on the obliquity angle (e.g. quasi-parallel models), 
the modulation of the shock modification with the obliq- 
uity angle is not appreciable in the presence of ejecta 
clumping (see Fig. |9l). In other words, our model pre- 
dict that the ejecta clumping can wash out the CR back- 
reaction effects on the separation between the forward 
shock and the contact discontinuity. On the other hand, 
our simulations have shown that the effects of back- 
reaction of accelerated CRs can still be visible on the 
azimuthal profile of plasma density which shows local 
maxima where the acceleration of CRs is the largest (the 
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Fig. 8. — Azimuthal profiles of the ratio of the forward shock 
radius to the contact discontinuity radius Rf^ /Red for models with- 
out back-reaction of accelerated CRs and with ejecta clumping 
and initial ejecta density profile either exponential (left panels) or 
power- law (right panels). The figure shows the profiles derived 
from models with maximum density perturbation i^max = 1.5 (up- 
per panels), I'max = 2.5 (middle), and i^max = 5.0 (lower), and with 
initial size of the clumps 1% (red lines) and Rj 2% (blue lines) 
of the initial diameter of the remnant DsnrO ■ The gree n line marks 
the p rofile derived from the observations of SN 1006 IIMiceli et al.l 
[2003) . 
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9. — As in Fig.|4]for a model accounting for the back-reaction 
of accelerated CRs only (run PL-QPAR-G1.3) and for models in- 
cluding both the CR, particle acceleration and the ejecta clumping 
(runs PL-C3.5-D1-QPAR-G1.3, and PL-C3.5-D2-QPAR-G1.3). In 
all the cases, the minimum ■y^g = 4/3. The green line marks 
the p rofile derived from the observations of SN 1006 ( Miceli ct sJ} 
\20m . The gray stripes mark the regions where the acceleration of 
CRs is the largest. 



plasma compressibility being the highest there). 
To make a more quantitative comparison between the 
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Fig. 10. — Median values of the ratio of the forward shock ra- 
dius to the contact discontinuity radius R{a/Rcd versus the max- 
imum density perturbation i^max for models either with (red and 
blue symbols for initial size of the clumps 1% and 2% of the 
remnant diameter DsnrOi respectively) or without (black symbols) 
ejecta clumping, and for models including the back-reaction of ac- 
celerated CRs (empty and crossed symbols for the quasi-parallel 
and isotropic cases, respectivel y) . The grey region marks the range 
of values observed in SN 1006 IIMiceli et al.ll2009h . 

model results and the observations, we derived the me- 
dian values of Ris/Rcd for each of the models in Table [1] 
and for the observed profile. Fig. [10] shows the median 
values of Rfs/Rcd versus the maximum density pertur- 
bation z/max for models accounting for only one of the 
effects considered in this paper (either back-reaction of 
accelerated CRs or ejecta clumping) and for models in- 
cluding both physical effects. We found that: the larger 
the size of initial clumps of ejecta, the lower the value of 
the median ratio; the higher the initial density perturba- 
tion, the lower the value of the median ratio. The back- 
reaction of accelerated CRs slightly reduces the value of 
the ratio in models accounting for the clumpy structure 
of the ejecta (empty symbols in Fig. [TU]) unless the en- 
ergy losses to CRs are large with an effective adiabatic 
index w 1.1 and ubiquitous at the forward shock (as in 
the isotropic injection, see models EX-ISO-Gl.l and PL- 
ISO-Gl.l; crossed symbols in Fig. [T0|. 

4. SUMMARY AND CONCLUSIONS 

We investigated the role of ejecta clumping and back- 
reaction of accelerated CRs on the evolution and mor- 
phology of young Type la SNRs and, in particular, on 
determining the observed separation between the for- 
ward shock and the contact discontinuity and the high 
occurrence of protrusions. To this end, we developed a 
3D MHD model describing the expansion of the remnant 
through a medium with nonuniform interstellar magnetic 
field, including consistently the back-reaction of acceler- 
ated CRs and the initial clumpy structure of the ejecta. 
We explored two complementary cases in which one or 
the other of these physical processes is turned either on 
or off in order to identify its effects on the remnant evo- 
lution and morphology. Then we compare d the model 
resul ts with the observations of SN1006 (|Miceli et al.l 
|2009() . Particular attention has been devoted to perform 
simulations with sufficient spatial resolution to capture 
the details of the evolution of the clumps of ejecta, ex- 
ploiting the adaptive mesh refinement capabilities of the 
FLASH code. 



As expected, we found that the acceleration of CR par- 
ticles makes the shell of shocked ISM thinner at the for- 
ward shock, thus reducing the separation between the 
forward shock and the contact discontinuity. Any de- 
pendence of the back-reaction of accelerated CRs on the 
obliquity angle should be evident as a modulation of the 
azimuthal profile of the ratio of the forward shock ra- 
dius to the contact discontinuity radius Rfs/Rcd- In the 
case of SN 1006, the comparison of the modelled pro- 
files with those observed shows that the back-reaction 
of accelerated CRs may reproduce the observations only 
if the energy losses to CRs are extreme (i.e. the effec- 
tive adiabatic index is joS ~ 1-1) and independent of 
the obliquity angle (i.e. the effects of CR acceleration 
are ubiquitous at the forward shock). In addition, the 
simulations have shown that the large compression ratio 
due to the acceleration of CR particles has no signif- 
icant effect on the growth of R T instabilities, in agree- 
ment with previous st udies (e.g. iBlondin fc Ellisonll2001l : 
iFraschetti et al.ll2010l : IWang||2011[ ). As a resuff, the rem- 
nant outline is only slightly perturbed by the instabilities 
with very few (if any) occurrence of protrusions near the 
epoch of young SNRs as SN 1006 or Tycho's SN R even 
with very efficient acceleration of CRs (see also IWangj 
1201 It) . This fact contrasts with the evidence of several 
protr usions observed in SN10 06 (e.g. [RakojTOkE^^ 
20Tll) and Tvcho's SN R, (e.g. iHwang fc Gotthell [l997l : 



Velazquez etallll998t ). 



On the other hand, the clumpy structure of the ejecta 
can have important consequences on the structure of the 
intershock RT mixing and on the final morphology of the 
remnant. In particular, we found that the ejecta clumps 
with the higher density contrasts approaching the con- 
tact discontinuity enhance the growth of RT instabili- 
ties; RT fingers can easily reach the forward shock and 
ejecta clumps can be found very close to, or even beyond, 
the average shock radius with no need to invoke any CR 
back-reaction at all to explain this phenomenon. As a re- 
sult, the separation between the forward shock and the 
contact discontinuity can be significantly reduced, de- 
pending on the size and density contrast of the clumps. 
In particular, we found that the larger the size of initial 
clumps of ejecta and/or the higher their density contrast, 
the shorter the width of the interaction region between 
the forward shock and the contact discontinuity. The 
modelled azimuthal profile of i?fs/i?cd is fairly uniform 
as observed in SN 1006; the comparison of the model re- 
sults with the observations of SN 1006 showed that the 
observed profile of Rfs/Rcd can be reproduced by mod- 
els with a maximum density perturbation of ejecta fc'max 
ranging between 2.5 and 5, and with initial size of ejecta 
clumps of the order of 2% of the initial diameter of the 
remnant. We also found that the remnant outline can 
be significantly perturbed by the enhanced RT fingers 
and, in case of high density contrasts and large size of 
the clumps, several protrusions can characterize the mor- 
phology of the remnant at the age of SN 1006. Our study 
supports the idea that enhanced RT mixing due to ejecta 
clumping can be responsible for the filamentary struc- 
tures and bumps seen on the outlines of young SNRs as 
SN 1006 and Tycho's SNR. 

Finally, our analysis has shown that the ejecta clump- 
ing, if present, may wash out the effects of back-reaction 
of accelerated CRs on the separation between the for- 
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ward shock and the contact discontinuity. In particular, 
if the CR acceleration efficiency depends on the obliquity 
angle as, for instance, in the quasi-parallel scenario, the 
modulation of the shock modification with the obliquity 
angle may be not appreciable in the presence of ejecta 
clumping. We conclude therefore that, in general, the 
separation between the forward shock and the contact 
discontinuity is not a reliable diagnostic tool for study- 
ing the CR shock modification. 

On the contrary, our model predicts that the effects of 
back-reaction of accelerated CRs can still be apprecia- 
ble on the azimuthal profile of plasma density. In fact 
our simulations have shown that, even in the presence 
of ejecta clumping, the density profile has local maxima 
where the acceleration of CRs is the largest (the plasma 
compressibility being the highest there). Also, due to 
the enhanced plasma compressibility, the magnetic field 
strength can reach values of ss 50 — 70 /iG where the 
CR acceleration is the largest (see Sect. 13. ip . that are 
comparable wi th those obse rved in the X-ray rims of 
SN1006 (e.g. iBerezhko et al 2003,, 120091: lAcero et al.l 
[Mot iPetruk et al.ll201lLl20Ti) . It is interesting to note 
however that similar values of magnetic field strength can 

http: // www. astropa. unipa. it/ progetti_ricerca /HPC/ index. html 



also be reached locally in ejecta clumps close to the for- 
ward shock with no need to invoke any CR back-reaction, 
but as a result of the propagation of the clumps through 
the intershock region (see Fig. [5]). 
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Fig. Al. — As in Fig.[TO]for models starting as early as either 2 yr (empty symbols; runs EX-C1.5-D1-2YR and EX-C5.0-D2-2YR) or 
~ 10 yr (filled symbols; runs EX-C1.5-D1, EX-C5.0-D2) after the SN explosion. Red and blues symbols show the results for models with 
initial size of the clumps 1% and ^ 2% of the remnant diameter DsnrO > respectively. The grey region marks the range of values observed 
in SN 1006 lIMiceli et al.|[200l l. 

APPENDIX 

DEPENDENCE OF THE RESULTS ON THE INITIAL CONDITIONS 

We checked the dependence of the results on the initial conditions and in particular on the time when the clumpy 
structure of the ejecta is initialized. To this end, we performed two additional simulations (runs EX-C1.5-D1-2YR and 
EX-C5.0-D2-2YR in Tabled]) with the same parameters of runs EX-C1.5-D1 and EX-C5.0-D2 but starting as early as 
« 2 yr after the SN explosion (i.e. the initial spherical remnant has radius Rfso = 0.125 pc) instead of « 10 yr (with 
Rfso — 0.5 pc). In other words, we checked the dependence of the results on the initial conditions for the limit cases 
considered in this paper, namely models with a clumpy structure of the ejecta characterized either by clumps with 
small size and low density perturbations (run EX-C1.5-D1) or by clumps with large size and high density perturbations 
(run EX-C5.0-D2). For'runs EX-C1.5-D1-2YR and EX-C5.0-D2-2YR, we used 13 nested levels of resolution in the 
automatic mesh derefinement scheme to keep the same spatial resolution as the other simulations discussed here (i.e. 
172 zones per initial radius of the remnant); in this case the effective mesh size was 32768 x 32768 x 32768. 

From the additional simulations, we derived the azimuthal profile of the ratio of the forward shock radius to the 
contact discontinuity radius Rfs/Rcd as done for the other runs (see Fig. [8|); then, from these profiles, we derived 
the median values of i?fs/i?cd and compared them with those derived from runs EX-C1.5-D1 and EX-C5.0-D2 (see 
Fig. lAip . In both cases analyzed, with ejecta clumps with either small size and low density perturbations or large size 
and high density perturbations (namely the two limit cases explored in this paper), we found that the median values 
of Rfs/Rcd derived from models with different initial ages are consistent within the error bars, the value being slightly 
lower (higher) in the model with initial age of the remnant tgnro = 2 yr than in the model with ismo = 20 yr when 
the clump size is 1% (2%) and the maximum density perturbation is = 1.5 (^ = 5). We conclude therefore that 
the results presented here do not depend on the initial age of the simulated remnant. Indeed our results undoubtedly 
show that the average separation between the contact discontinuity and the forward shock strongly depends on the 
clumpy structure of the ejecta and, in particular, on the size and density contrasts of the clumps. 



